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Abstract: Increasing demand for natural rubber (NR) calls for an increase in latex yield 
and also an extension of rubber plantations in marginal zones. Both harvesting and abiotic 
stresses lead to tapping panel dryness through the production of reactive oxygen species. 
Many microRNAs regulated during abiotic stress modulate growth and development. The 
objective of this paper was to study the regulation of microRNAs in response to different 
types of abiotic stress and hormone treatments in Hevea. Regulation of MIR genes differs 
depending on the tissue and abiotic stress applied. A negative co-regulation between 
HbMIR398b with its chloroplastic HbCuZnSOD target messenger is observed in response 
to salinity. The involvement of MIR gene regulation during latex harvesting and tapping 
panel dryness (TPD) occurrence is further discussed. 
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1. Introduction 

Hevea brasiliensis is the only source of natural rubber (NR), which is produced in latex cells in the 
phloem. Latex harvesting consists in tapping the soft bark of rubber trees. Ethephon is applied to 
stimulate production. Increasing demand for NR calls for an increase in latex yield, but also an 
extension of crops to marginal zones. Such marginal zones are subject to stronger and more 
detrimental abiotic constraints for latex production (positive cold, frost-prone, drought, etc.). Exposure 
to abiotic stress in addition to latex harvesting stress affects latex production and tree productivity [1]. 
Laticifer cells are the site of numerous types of stress that lead to oxidative stress. The plant's first 
response to wounding is to produce jasmonate [2]. The wound and associated jasmonate production 
also cause oxidative stress in the plant and activate antioxidant systems to overcome it [3]. Oxidative 
stress, whether produced in response to tapping and/or environmental constraints, may lie behind 
problems linked to latex flow [4]. It is reflected in an imbalance between ROS (reactive oxygen 
species) and the detoxification systems existing in laticifers. ROS then attack the membranes of cell 
organelles, causing the release of coagulant factors existing in the lutoids, which gives rise to in situ 
coagulation of rubber particles. This complex physiological disorder is called TPD (Tapping Panel 
Dryness). This phenomenon remains reversible up to a certain limit. It amounts to laticifer fatigue, 
called dry cut. In extreme cases, halted flow is followed by the degenerescence and death of the 
laticifers and, consequently, severe cortical necrosis (Brown Bast) [5]. Annual losses due to TPD are 
estimated at 10% to 40% [1]. 

MicroRNAs (miRNAs) form a small class of non-coding RNAs of 21 to 24 nucleotides involved in 
numerous physiological processes in plants including response to biotic and abiotic stresses through 
post-transcriptional regulation. miRNA biogenesis requires several stages and several enzymes from 
MIR gene transcription to the generation of mature miRNA [6-10]. The mature miRNA enters the 
complex called RISC (RNA-induced silencing complex) thereby enabling regulation of the target gene. 
In the RISC complex, the mature miRNA associates with proteins of the Argonaute family (AGO) to 
direct the regulation of the target gene [11,12]. In the RISC complex, the mature miRNA sequence 
links up with the sequence of the messenger RNA through perfect or imperfect sequence homology. 
This action causes a silencing phenomenon in plants, called PTGS (post-transcriptional gene silencing) [13]. 
In plants, most target messenger RNAs contain a single site complementary to the miRNA, 
corresponding to the cleavage site [9], suggesting that plant microRNAs use RNA cleavage rather than 
preventing translation. Long microRNAs are taken in charge by Dicer-like 3 (DCL3), associate with 
AG04 and then guide methylation of the DNA of target genes in rice and mosses [14,15]. 

Many microRNAs regulated by abiotic stress have been identified in the model species (for review, 
see [16]). The first involvements of microRNAs in response to stress were described by Jones Rhoades 
and Bartel who, in Arabidopsis thaliana, predicted genes targeted by microRNAs such as superoxide 
dismutase, the laccases and ATP sulfurylases (APS) [17]. MicroRNA expression was then studied in 
different contexts of abiotic stress in numerous species (for review, see [16]). Most microRNAs 
conserved between species target transcription factors, such as miR/target gene pairs miR156/SBP, 
miR159/319/MYB-TCP, miR160/ARF, miRI66/HDZIPIII, miRl 69/NFY subunit 

MicroRNAs are also involved in the regulation of antioxidant activities and particularly in the 
regulation of reactive oxygen species-scavenging enzymes. The production of ROS is linked to abiotic 
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and biotic stress status in plants [18-20] and the excess ROS can lead to several types of cell damage, 
such as peroxidation of membrane lipid compounds, degradation of polysaccharides, denaturing of 
enzymes and lesions in DNA [21]. The miR398 targets the CSD1 and CSD2 genes encoding cytosolic 
and chloroplastic CuZnSOD in Arabidopsis thaliana [22]. The level of expression of these microRNAs 
and of their target is regulated in response to stress, involving a modulation of growth and 
development during stress. For example, miR398 is under-expressed in response to stress which leads 
to an accumulation of CuZnSOD enabling detoxification of an excess of reactive oxygen species in 
Arabidopsis thaliana [22]. 

To date, 28 sequences of MIR are accessible in miRBase [23-25] for the rubber tree. They represent 
10 conserved families and 14 newly identified families potentially 7/evea-specific miR [26,27], added to 
which four additional precursors have been identified recently [28]. Our previous studies on mature trees 
demonstrated the functionality of miRNA biogenesis in the latex cells. Moreover, we highlighted the 
reduction of small RNA size from 24 to 21 nucleotides in TPD-affected trees, which could not be 
explained by a general RNA degradation [28]. From this previous work, we hypothesized that the 
change in size could also be effective for mature miRNAs and consequently involved to independent 
miRNA pathway, one leading to targeted transcript cleavage and the other to methylation of targeted 
genes (for review [29]). Among all pre -miRNAs tested, only the relative Hbpre-MIR159b abundance 
was up-regulated between TPD-affected and healthy trees. Targets of miR159 were predicted and 
might be involved in various activities such as rubber biosynthesis, antioxidant activity and 
transcription regulation activity [26,28,30]. 

The objective of this paper was to study the regulation of microRNAs in response to different types 
of abiotic stress (cold, saline stress, wounding) or hormone treatments (ethylene and methyl jasmonate). 
Given some mature miRNAs could be produced by several MIR genes, this study focused on the 
expression analysis pattern using real-time RT-PCR with specific primers for each MIR genes. The 
relative accumulation of 19 premiRs was analyzed by real-time RT-PCR on various tissues from 
juvenile plant material in order to facilitate the application of treatments. Negative co-regulation of 
Hbpre-miR398b and its target chloroplastic CuZnSOD is revealed in response to saline stress. 

2. Results and Discussion 

2.1. Relative Accumulation of premiR in Response to Abiotic Stress in in Vitro Plantlets 

Relative accumulation of premiR of 20 selected MIR genes was analysed in different tissues 
(leaf, bark and root; Table 1). In control plants, the conserved MIR genes were generally very weakly 
expressed such as HbMIR156, HbMIR319 and HbMIR398a and HbMIR398b (from 6.78 x 10" 4 in roots 
for HbMIR398a to 9.22 x 10" 2 for HbMIR319 in bark). By contrast, the HbMIR408a and HbMIR408b 
genes were strongly expressed (from 8.47 x 10 to 3.54). Among the MIR genes newly identified in 
Hevea brasiliensis, HbMIR6482, HbMIR6483 and HbMIR6485 were also strongly expressed 
(from 2.20 for HbMIR6482 to 7.44 x 10 for HbMIR6485 in leaves). 

In the same in vitro plantlets, statistical analyses revealed a difference in gene expression at 
tissue level (Table 1). Indeed, three genes HbMIR166b, HbMIR319 and HbMIR6482 had a relative 
differential expression level in leaves, bark and roots under untreated conditions. HbMIR166b and 
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HbMIR319 were more strongly expressed in bark while HbMIR6482 was more strongly expressed in 
roots. For these three genes, the lowest expression level was found in leaves (Table 1). In addition, 
seven genes significantly displayed a relative differential expression level, namely HbMIR159a, 
HbMIR398a, HbMIR398b, HbMIR408b, HbMIR6483, HbMIR6484 and HbMIR6485. The MIR398c 
transcript was only present in trace form under the majority of conditions. 

Table 1. Relative premiR abundance of 13 conserved MIR genes and 7 MIR genes newly 
identified in Hevea brasiliensis in the leaves, bark and roots of one-year-old in vitro 
plantlets in response to cold and saline (NaCl) treatments compared with non-treated plants 
(C: control), nc: not calculated. Data with the same letter are not significantly different at 
the level of 5%. 
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Table 1. Cont. 
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The regulation of MIR genes in response to cold and salinity stresses was shown as the ratio 
between the mean of three treated in vitro plantlets and three untreated in vitro plantlets (Table 1). 
An analysis of the relative expression of the MIR genes showed that 8 MIR genes were significantly 
regulated in response to cold (Table 1). In leaves, both HbMIR159 genes displayed a significant 
increase in the relative number of transcripts in response to cold while the relative expression of 
HbMIR6482 was greatly inhibited. In bark, the relative expression of 3 MIR genes (HbMIR166b, 
HbMIR476 and HbMIR6482) was repressed while that of HbMIR159b was slightly stimulated. 
In roots, 5 MIR genes (HbMIR159a, HbMIR408a, HbMIR408b, HbMIR6482 and HbMIR6484) 
displayed a significant reduction in transcripts in response to cold. The drop in temperature led to a 
drop in transcripts for HbmiR6482 in the 3 tissues. HbMIR159a displayed an expression profile that 
was antagonistic between leaves and roots. With the exception of HbMIR159a and HbMIR159b, the 
significantly regulated MIR genes were repressed in response to cold stress. 

The in vitro plantlets were also subjected to saline stress by watering with a 300 mM NaCl solution 
rather than water twice a day for 24 h. In response to NaCl, 16 of the 20 genes displayed a significant 
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repression in their level of relative transcript abundance. Only HbmiRnl2 was highly induced (4 to 6 times) 
in most tissues. 

2.2. Relative Accumulation of preMIR in Response to Ethylene (ET), MethylJasmonate (MeJA) and 
Wounding in Three-Month-Old Epicormic Shoots 

Based on previous studies on budded plant material [31], the effect of ET, MeJA and wounding was 
studied on three-month-old epicormic shoots. In control plants, the relative level of premiR was highly 
accumulated for HbMIR6485 and HbMIRnB (1.60 x 10 and 3.92 x 10 2 respectively) and at low level 
for HbMIR156, HbMIR319 and HbMIR398a {AM x 10" 4 , 1.15 x io" 2 and 3.61 x 10" 4 , respectively; 
Table 2). This result was similar to the basal expression profile found in the in vitro plantlets. 
7 MIR genes had a relative expression level that differed significantly between bark and leaves (Table 2). 
They were HbMIR159a, HbMIR159b, HbMIR166b, HbMIR319, HbMIR396, HbMIR398a and 
HbMIR6485. In comparison, the aforementioned genes were also regulated differently depending on 
the tissues in the in vitro plantlets, apart from HbMIR396. Again, the MIR398c transcript was only 
present in trace form under the majority of conditions. 

Table 2. Relative premiR abundance of 13 conserved MIR genes and 7 MIR genes newly 
identified in Hevea brasiliensis in the leaves and bark collected on 3 -month-old epicormic 
shoot of budded plants in response to ethylene (ET), methyl jasmonate (MeJA) and 
wounding (W) compared with non-treated plants (C: control), nc: not calculated. Data with 
the same letter are not significantly different at the level of 5%. 
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3.52 
0.56 


0.08 
0.24 


1.35 
0.60 


0.98 
0.39 


HbpremiRnl 2 


leaves 
bark 


6.19 x I0~ 2a 
6.13 x I0" 2a 


1.23 
1.20 


0.65 
0.58 


1.51 
0.84 


0.12 
0.96 


7.12 
1.42 


0.13 
0.44 


HbpremiRnl 3 


leaves 
bark 


3.16 x 10 2a 
3.92 x I0 2a 


1.23 
1.85 


0.39 
0.16 


1.24 
0.91 


0.55 
0.55 


1.41 
1.85 


0.96 
0.63 



The regulation of MIR genes in response to ethylene, MeJA and wounding was shown as the ratio 
between the mean of three treated in vitro plantlets and three untreated in vitro plantlets. 3 MIR genes 
were significantly regulated by the ethylene treatment (Table 2). In leaves, only the relative expression 
of the HbMIR398a gene was repressed. In bark, 2 MIR genes were negatively regulated by ethylene, 
and the expression of HbMIR319 and HbMIR396 was inhibited (Table 2). Application of 
methyl jasmonate induced regulation of the relative expression of 6 MIR genes in leaves and not in 
bark (Table 2). The HbMIR159a, HbMIR159b, HbMIR166a, HbMIR476 and HbMIR6485 genes 
displayed a relative increase in transcripts, while HbMIR408b displayed a significant large drop in 
transcripts in response to external methyl jasmonate (Table 2). In response to wounding, 5 genes were 
significantly regulated (Table 2). In leaves, HbMIR408a and HbMIR408b displayed a large reduction 
in their relative expression levels while HbMIR6482 significantly showed an increase in relative 
abundance for these transcripts (Table 2). In bark, the HbMIR166b and HbMIR319 genes displayed a 
significant drop in their relative transcript abundance (Table 2). It should be noted that the HbMIR319 gene 
was regulated in the same way in response to both ethylene and wounding. 

The response to jasmonate in leaves seemed to differ from the other abiotic stresses. In fact, the 
expression of 5 MIR genes (MIR159a, MIR159b, MIR166, MIR476 and MIR6485) was stimulated 
during the treatment. MiR159 and miR476 are known to target type MYB transcription factors [17,32] 
and PPRs (pentatricopeptide repeat proteins) [33], respectively. These targets were not identified 
during the prediction of targets of miR159 and miR476 in Hevea. In contrast, miR166 is known to 
target proteins of the HD-ZIPIII family in Arabidopsis and also in Hevea [17,26]. 
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The expression of MIR genes was differentially regulated depending on the stress and the tissue 
(bark and leaf), suggesting specific response pathways for each type of stress (Figure la,b). 

Figure 1. Venn diagram representing the regulation of the MIR genes common to all the 
treatments tested, (a): in leaves; and (b): in bark. 




BARK 



(b) 
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In maize, conversely, the comparison of expression pattern of mature miRNAs in plant subjected to 
saline stress and a viral infection lead the authors to suggest some cross-talk between abiotic and biotic 
stress [34]. In addition to differential expression depending on the tissue and stress applied, some 
differences between species are worth noting. In Hevea brasiliensis, the response to cold stress differs 
from that found in Arabidopsis thaliana [35]. In Arabidopsis, miR166b is accumulated in response to 
cold stress whereas in Hevea, repression of the HbMIR166b gene is found. Inhibition of the relative 
expression of the HbMIR476 gene in Hevea bark shows an expression profile similar to that in poplar [36]. 
For saline stress, regulation of the expression of microRNAs appears to be specific to each species. 
In Hevea, saline stress led to global repression of the expression of all the MIR genes studied, 
except HbMIRnl2. In this case, the result suggests that response mechanisms were brought into play 
by blocking the repression of genes targeted by microRNAs. By contrast, in Arabidopsis, accumulation 
of the mature microRNAs miR156, miR159, miR319 and miR396 and a reduction in the quantity of 
miR398 are seen [37]. The quantity of miR396 and miR156 is reduced in rice and maize 
respectively [38]. In poplar, the quantity of miR530a, miR1445, miR1446a-e, miR1447 and miR1711 
is reduced while that of miR482.2 and miR1450 is increased [36]. 

The difference observed between Hevea and other plants could be explained by post-transcriptional 
regulations acting between pri-MIR transcript and the release of mature miRNA [39]. Moreover, 
most of studies reported on total mature miRNAs using Northern-blot hybridization highlighting 
the final products of several MIR genes. In our study, the use of real-time RT-PCR allowed 
identifying the specific expression of each MIR gene isoform. The microRNAs newly identified in 
Hevea brasiliensis have a high level of expression compared to that of conserved microRNAs under 
normal conditions. Different elements tend to show that the new microRNAs (young miRNAs) have 
little or no function [40]. Although conserved microRNAs are generally more abundant than the new 
microRNAs [41], an abundance of the latter does not mean greater regulation of the targets. Regulation 
of the targets of new microRNAs is not affected in mutant transgenic plants in the miRNA biogenesis 
pathway. The new microRNAs would seem not to be integrated yet in the regulation networks or may 
only function in a precise spatio-temporal context, or during response to a specific type of biotic or 
abiotic stress [40]. 

2.3. Analysis of the Co-Regulated Expression of the HbMIR398a, HbMIR398b anJHbMIR398c 
Genes and Their Putative Target Gene Chloroplastic HbCuZnSOD in Response to Abiotic Stress and 
to Hormone Treatments 

Three miR/target pairs have been experimentally validated in the rubber tree (miR156/Squamosa 
promoter binding protein, miR160/ARF and miR398/chloroplastic CuZnSOD) [26]. We chose to 
monitor expression of the MIR398/ chloroplastic CuZnSOD pair. HbmiR398 is generated by three 
HbMIR398 genes (3 genes, a, b and c; Figure 2). 

In order to monitor expression of the chloroplastic HbCuZnSOD target, two primer pairs were 
designed, one flanking the cleavage site and the other at 3' UnTranslated Region (3TJTR) (Table 3). 
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Figure 2. Stem-loop structure of three precursors of MIR398 in Hevea brasiliensis . 
The mature microRNA sequence is shown in dark. 
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Table 3. List of primers and their sequences used for MIR gene expression analysis by 
RT-PCR. All primers are presented from 5' to 3' end. 



Pre-microRNA 


Forward primer 


Reverse primer 


PCR 
efficiency 


Hbpre-miR156 


TGGTGATGTTGTTGACAGAAG 
ATAGAGAGC 


GCACAAAGGAGTGAGATGCA 
GAGTCC 


1.79 


Hbpre-miRl 59a 


GGTTAAGAAGTGGAGCTCCTT 
GAAGTC 


GCTCCCTTCAATCCAAACAAG 
GATC 


1.958 


Hbpre-miR159b 


GTGGAGCTCCTTGAAGTCCAA 
TAGAGG 


AGAGCTCCCTTCAATCCAAAC 
AAGG 


1.881 


Hbpre-miRl 66a 


TTCTTTTTGAGGGGAATGTTG 
TCTGG 


GGAATGAAGCCTGGTCCGAG 
GAG 


1.820 


Hbpre-miRl 66b 


GGGGAATGTTGTCTGGTTCGA 
TG 


TCAAATCAAACCCTGTTGGGG 
G 


1.738 


Hbpre-miR319 


CCAGTCACGGTGGGCAATGGG 


GGAGCTCCCTTCAGTCCAAGT 
ACAGG 


1.847 


Hbpre-miR396 


TGACCCTCTTCGTATTCTTCCA 
CAGC 


CCCACAGCTTTATTGAACCGC 
AAC 


1.782 


Hbpre-mir398a 


TGAGAACACAGGTGTTTTGGC 
TACC 


GTGCTCCAAAGGGGTGACCTG 
AG 


1.879 


Hbpre-mir398b 


ACCTGAGATCACATGTGGACA 
CCC 


GCGGTGGAGGAGAGCCCAG 


1.939 


Hbpre-mir398c 


TGGCCACCCTCACATGTTCCC 


CCGGCAGGGGTGACCTGAG 


1.965 


Hbpre-miR408a 


ACTGGGAACAGGCAGAGCAT 
GG 


GCCACAAGCCAGGGAAGAGG 
C 


1.723 
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Table 3. Cont. 



Pre-microRNA 


Forward primer 


Reverse primer 


PCR 
etliciency 


Hbpre-miR408b 


GACATACAAAGACTGGGAAC 
AGGCAG 


GCCACAAGCCAGGGAAGAG 
GC 


1.792 


Hbpre-miR476 


GCCTTGTATGTTTCATTTAGTA 
A 1 CC 1 1 C 1 


GATAATCCTTCTATGCAAAGT 

/~i r T"~T"~T"~T' A T/° Z - ^ 

CI 1 1 1A1GC 


1.732 


Hbpre-miR6482 


ACCAGGAACTGGTATCAACCC 
AGC 


TGCTACCAATGAATCGGACCC 

A Z~^ /~1 

ACC 


1.837 


Hbpre-miR6483 


CAGTAAATAGCAGTATCGTGG 

ATA 

A I AGGG 


GTCCAATCATTGATCCTGAAA 

a TT ,r rr ,r r a z~^ 

All 1C1AC 


1.828 


Hbpre-miR6484 


TGGATTGGAGCCCAATACTGT 

Z~^ A (~~* 


CTGCTCCATTGATTTTACCATC 
1 Al CjL 


1.873 


Hbpre-miR6485 


A r^C^T 1 A A A A A A Z" 1 Z" 1 A 

ALL 1 AGGA 1 G 1 AGAAGAGCA 
TAAC 


A /"^T A (~~* A T^Z - ^ A Z - ^ T -1 Z - ^ Z - ^ A ~F A T"' A ~F A Z~^ 

AC 1 ACA 1 GAG 1 GGA 1 A 1 A 1 AG 
GAATCC 


1.787 


Hbpre-miRnl 1 


G 1 A 1 CAACGCAGA 1 G 1 GCCGC 
C 


CCCCAGCCAAACTCCCCACC 


1.828 


Hbpre-miRnl 2 


AGC 111 CACCCAA 1 AACC 111 
LiUALr 1 


Z - ^ Z - ^ - FZ~~' r ~^' - FZ _, Z - ^ A A TTTT'T A TPP A A 

GC 1 C 1 1 CCAA 1 1 CC 1 A 1 CCAA 
a r,Tr,r,T 

PiXJ 1 UU 1 


1.78 


Hbpre-miRnl 3 


TGTGTTGGCCTTCGGGATCGG 


CGAATGCCCCCGACTGTCCC 


1.889 


Hb-RH2b 


GAGGTGGATTGGCTAACTGAG 
AAG 


GTTGAACATCAAGTCCCCGAG 
C 


1.68 


HbCuZnSOD 
(flanking miRNA site) 


GCTCTATCTCTCGCCGCCGCC 
TCC 


CCGCAATTGTTGCTTCTGCC 


1.785 


HbCuZnSOD (3'UTR) 


TGGCAGAAGCAACAATTGCGG 


GCAGGGAACAATGGCTGCC 


2 



In response to abiotic stress, the transcripts of the chloroplastic HbCuZnSOD were accumulated 
significantly in roots in response to cold, in bark and roots in response to saline stress, and particularly 
after 24 h of treatment revealing that these transcripts are not cleaved by miR398 (Table 4, Tables SI 
and S2). When the relative abundance of chloroplastic HbCuZnSOD transcripts was visualized 
(cleaved and non-cleaved with the primers at 3'UTR), the relative expression level for 
chloroplastic HbCuZnSOD transcripts decreased in leaves after 24 h of saline stress. However, that 
relative expression level rose significantly in response to saline stress in roots after 24 h of saline 
stress, and in wounded leaves (Table 4). 

The relative accumulation of premiRs was analyzed for the three MIR398 genes. The MIR398c 
transcript was only present in trace form under the majority of conditions. The relative expression of 
the MIR398a and MIR398b genes was significantly repressed in response to stress. In fact, in leaves, 
the relative expression of the MIR398a gene was repressed in response to 24 h of saline and ethylene 
treatment. Likewise, the relative expression of the MIR398b gene was repressed in bark and roots after 
24 h of saline stress treatment. 

When the relative expression level for the HbMIR398 gene was compared with the expression data 
for the target, the significant reduction in the expression of the HbMIR398b gene in bark and roots in 
response to 24 h of saline stress was accompanied by a significant increase in the relative expression 
level for the non-cleaved chloroplastic HbCuZnSOD transcripts (Figure 3). 
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Table 4. Expression profiles of the HbMIR398a, HbMIR398b and HbMIR398c genes and 
their chloroplastic HbCuZnSOD target by real time-PCR in leaves, bark and roots of 
one-year-old Hevea brasiliensis in vitro plantlets in response to a positive cold, NaCl, 
ethylene, methyl jasmonate (Me J A) and wounding treatments. The statistically significant 
ratios (p-value < 0.05; Tables SI and S2) with a value over 1 (in red) and under 1 (in 
green), represent over-expression and under-expression, respectively, in response to the 
treatments. The non-significantly regulated genes are shown in yellow (nc: not calculated). 



Ratio of relative transcript accumulation (Treated/Control) 



Treatment 


Tissue 


chloro CuZnSOD chloro CuZnSOD 
(flanking miRNA site) (3'UTR) 


preMIR398a 


preMIR398b 


preMIR398c 




leaf 


2.71 


3.76 


0.34 


0.95 


nc 


Cold 


bark 


8.98 


8.90 


0.23 


0.52 


nc 




root 




1.46 


0.47 


0.62 


nc 




leaf 


0.84 


0.29 


0.21 


0.87 


nc 


NaCl 


bark 




5.21 


0.49 




nc 




root 


37.45 




0.1 


0.15 


nc 


Ethylene 


leaf 


nc 


0.42 


0.22 


0.52 


0.91 


bark 


0.49 


0.32 


3.59 


15.19 


nc 


MeJA 


leaf 


1.47 


nc 


0.69 


0.09 


nc 


bark 


0.79 


1.00 


3.38 


0.61 


nc 


Wounding 


leaf 


2.61 




1.85 


0.05 


0.04 


bark 


0.85 


0.88 


2.99 


4.66 


nc 



Figure 3. Graph showing the expression values for the chloroplastic HbCuZnSOD and 
HbMIR398b genes in response to saline stress in bark and roots after 24 h of treatment. The 
expression values represent the mean and standard deviation of three treated or untreated 
plants, for each tissue. Data with the same letter are not significantly different at the level 
of 5%. 
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Co-expression of MIR398 and of the gene encoding chloroplastic CuZnSOD was analyzed in 
response to several treatments. This enzyme takes part in ROS detoxification through the transformation 
of two superoxide ions H2O2 and H2O. In Arabidopsis, miR398 targets cytosolic and chloroplastic 
superoxide dismutase (CSD1 and CSD2, [22]). Under oxidative stress conditions (heavy metal, strong 
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sunlight and methyl viologen treatments), negative regulation of CSD1 and CSD2 by miR398 is halted 
in leaves, leading to an accumulation of enzymes and a drop in toxic free radicals for the cell. 
Tolerance of oxidative stress is consequently increased [22]. Interestingly, the chloroplastic CuZnSOD 
is the only one to have been validated in Hevea for miR398, the cytosolic form was not predicted to be 
targeted by miR398 and the absence of cleavage was experimentally validated [26]. Some other targets 
have been predicted for this microRNA but have yet to be experimentally validated [26]. Co-regulation 
was only visible in response to saline stress. A reduction in the expression of the HbMIR398b gene 
was associated with an accumulation of chloroplastic HbCuZnSOD transcripts in bark and roots 
but not in leaves. It is noteworthy to specify that the bark of the young Hevea in vitro plantlets 
was also chlorophyllous. For root tissues, other plastid organelles might be present. In addition, the 
level of expression for the MIR398 genes was not always negatively correlated with the level of 
chloroplastic HbCuZnSOD transcripts, which has already been found in other species and suggests the 
existence of two regulation mechanisms, dependent on and independent of miR398 [42]. 

3. Experimental Section 

3.1. Plant Material and Treatments 

In vitro plantlets were regenerated by indirect secondary somatic embryogenesis from maternal 
inner integument of immature fruits from rubber clone PB 260 [43]. These in vitro plantlets were 
acclimatized and reared in the greenhouse for one year at a temperature of 28 °C. At that stage, the 
treatments were applied on three biological replicates. For cold stress, the plants were placed at 4 °C 
for 12 h. Saline stress was applied to the in vitro plantlets by watering with a 300 mM solution of NaCl 
rather than just water twice a day for the duration of the treatment (24 h). 

Three-month-old epicormic shoot from grafted plants were used for ethylene (ET), methyl-jasmonate 
(Me J A) and wounding treatments. Plants from the Hevea clone PB 260 were grafted onto rootstock 
GT1. This plant material included a mature growth unit corresponding to a three-month-old epicormic 
shoot after cutting back. The treatment was applied for 4 h on the whole plant at 8:00 a.m. 
For application of the ethylene and methyl jasmonate treatments, each plant was placed in a 
hermetically sealed plexiglas box with a volume of 300 L leaving the doors open for 24 h prior to 
treatment to limit the effects of transport or any wounds. For application of the ethylene treatment, the 
gas was injected to obtain an ethylene concentration of 1 ppm (300 uL of pure ethylene) in the box. 
For the methyl jasmonate treatment, 20 uL of methyl jasmonate solution (>95%) (Sigma, St. Louis, 
MO, USA) was diluted in 500 uL of absolute ethanol and placed on Whatman paper inside the box, 
in order to release the methyl jasmonate in gas form. For these two stress conditions, samples were 
taken 4 h after treatment. The plants used as controls for these treatments were also placed in a 
hermetically sealed box without gas injection. In that way, it was possible to eliminate the impact of 
the box when comparing the two treatments (control and treated). For the wounding treatment, the 
plants were wounded on their leaves by applying pressure every 2 cm with serrated tweezers and on 
the bark by making cuts with a scalpel 1 cm apart on the stem of the scion. 

Table 1 summarizes the different treatments applied to the rubber in vitro and grafted plants. Some 
leaf, bark and root samples were taken after 12 h for the in vitro plantlets, and only leaves and bark for 
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the budded plants whose root system was not clonal. All the samples were immersed directly in liquid 
nitrogen after being taken and were stored at -80 °C pending RNA extraction. 

3.2. Extraction and Purification of Total KN As 

The samples (leaf, bark and root) were ground in liquid nitrogen. The resulting fine powder was 
placed on ice in an extraction buffer (25 mM guanidium isothiocyanate 5M NaAc, 0.88% of sarcosine, 
0.9% of polyvinylpyrrolidone (PVP) and 1% of P-mercapto-ethanol). The samples were mixed for 30 s 
to optimize membrane lysis. An initial centrifugation (30 min at 4 °C at 10,000 g) was applied to 
precipitate cell debris. The supernatant was deposited on a caesium cushion (24 mM NaAc, 
45.6 mM CsCl). The tubes were then balanced and centrifuged for 20 h at 20 °C at 25,000 rpm in a 
Beckmann Coulter ultracentrifuge (L7, Beckmann Coulter, Brea, CA, USA) using the SW41 rotor 
(Beckmann Coulter). Only dense RNAs were able to cross the caesium chloride cushion. After 20 h, the 
supernatant was removed by pipetting, leaving 1 cm height of solution, which was eliminated by 
tipping up the tube. This precaution was taken to avoid contamination by genomic DNA. The pellet 
was dried for 5 min then rinsed in 70% ethanol then taken up in 200 uL of sterile water. The RNAs 
obtained were quantified then checked by electrophoresis (Mini-Sub Cell, Biorad, Richmond, CA, USA) 
on agarose gel (1% w/v, 1* Tris Acetate Ethylenediaminetetraacetic acid) and stained with ethidium 
bromide for visualization under ultraviolet. The RNAs were then aliquoted and stored at -80 °C. 

3.3. Checking for the Presence of Genomic DNA and DNAse Treatment 

A polymerase chain reaction (PCR) was carried out to check for the presence of genomic DNA in 
the samples of extracted total RNA. Amplification was carried out in a final volume of 50 uL. 
This was done with a volume of RNA corresponding to 50-100 ng to which was added a reaction 
mixture comprising 10x reaction buffer (10 mM Tris, pH = 8.3), MgCl 2 (1.5 mM), dNTP (0.04 mM), 
gene primers encoding actin (Forward: (5'-TCCATAATGAAGTGTGATGT-3', Reverse: 
(5'-GGACCTGACTCGTCATACTC-3', chosen astride an intron) (0.04 uM) and Taq polymerase 
(1 U/uL). A sample of Hevea genomic DNA was used as the positive control. The samples were 
incubated for 1 min and 30 s at 94 °C followed by 40 cycles, each composed of 30 s of denaturing at 
94 °C, 30 s of hybridization at 58 °C and 90 s of elongation at 72 °C, with the reaction ending on 
a final elongation of 10 min at 72 °C. The amplification products were deposited on an agarose gel 
(1% w/v, \x Tris Acetate Ethylenediaminetetraacetic acid) then stained with ethidium bromide. The 
RNA samples that proved to be contaminated by genomic DNA were treated by digestion with DNase 
(Turbo DNA-FREE kit from Ambion, Life technologies, Carlsbad, CA, USA). The reaction took place 
from 8 ng of total RNAs in a final volume of 40 uL composed of DNase buffer (10x, 4 uL, Ambion), 
DTT (0.1 M, 2 uL, Invitrogen, Life technologies, Carlsbad, CA, USA), RNase inhibitor (RiboLock, 
Fermentas, Burlington, ON, Canada, 40 U/uL) and TURBO DNase (2 U/uL). The samples were then 
incubated for 15 min at 37 °C. Next, a DNase inactivator (TURBO DNase Inactivation Reagent, 
Ambion) was added to the RNA samples. Enzyme inactivation took place for 2 min at ambient 
temperature by mixing the tubes. The supernatant containing the purified RNAs was recovered by 
centrifugation (12,000 g, 2 min, 4 °C). 



Int. J. Mol. Sci. 2013, 14 



19601 



3.4. cDNA Synthesis and Real-Time PCR 

A microgram of total RNAs was used to synthesize some cDNAs in a reaction volume of 20 uL 
using RevertAid™ M-MuLV reverse transcriptase (MBI, Fermentas, Burlington, Canada) following 
the supplier's recommendations. The primers used to amplify the stem-loop structures are listed in Table 2. 
The PCR conditions comprised a denaturing cycle at 95 °C for 2 min, then 45 cycles comprising 
denaturing at 95 °C for 20 s and a hybridization and elongation stage at 60 °C for 20 s. The PCR 
reaction took place in a volume of 6 uL containing 2 uL cDNA diluted 25 times, 1 uL of primers at 
5 uM, and 3 uL of 2x SYBR green mixture (LightCycler® 480 SYBR Green I Master, Roche Applied 
Sciences, Bale, Switzerland). All the primers were validated by generating a standard curve from the 
dilution in series of 8 points, in triplicate, to calculate the efficiency of each primer pair (Table 3). 
In addition, amplification specificity was checked by generating the melting curves that needed to 
display a unique peak, and all PCR products were checked by sequencing. Relative quantification was 
carried out with Light Cycler 480 software (version 1.5.0, Roche Applied Sciences, Bale, Switzerland), 
using RH2b [44] as the reference. 

3.5. Statistical Analysis 

The statistical analyses (analysis of variance followed by a Fisher comparison of means test) were 
done with XLSTAT (version 2011.4.02, Addinsoft, Paris, France) using expression data normalized 
by the LOG10(X) function. It was considered as up-regulation when the ratio was >1.0 and 
down-regulation when the ratio was <1.0. The p value corresponds to the Fisher test of the ANOVA. 

4. Conclusions 

Studying the regulation of MIR genes under harvesting stress (tapping, ethephon stimulation) 
and ROS-induced TPD in mature rubber trees remains difficult. However, some recent studies 
revealed that only one MIR genes is differentially regulated upon TPD. The MIR159b gene was shown 
up-regulated upon TPD occurrence [24]. The expression of this gene was increased in response to cold 
in leaves and bark, and in response to the jasmonic acid treatment in leaves of juvenile plantlets. 
Conversely, it was inhibited in all the tissues in response to salinity. 

In order to get a full understanding of mechanisms involved in latex production and TPD syndrome, 
a complete validation of miRNA/target messenger pairs is first needed using by high throughput 
"degradome" analysis [45]. Combination of analyses on juvenile and mature plant materials will help 
developing model of MIR gene regulations under abiotic stress and further characterization of the 
TPD-regulated miRNAs and their targets. 
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